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Abstract 

We present a QCD sum rule calculation of the strange-quark mass including 
four-loop QCD corrections to the correlator of scalar currents. We obtain 
rhsil GeV) = 205.5 ± 19.1 MeV. 
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1 .Introduction 



A precise determination of the values of the hght quark masses is of crucial practical 
importance for testing in an accurate way the predictions of the Standard Model. 
In particular, the knowledge of the strange quark mass is relevant for a better 
understanding of the low-energy phenomenology of QCD and for a precise prediction 
of the CP- violating parameter e'/ e in the framework of the Standard Model [|12|, 0, 



14. l". 



The ratios of the light quark masses can be determined in a model-independent 
way with the help of chiral-perturbation methods On the other hand, in order 
to obtain their absolute values, one has to resort either to the method of QCD sum 
rules or to lattice QCD ||16|, |30|. While no less fundamental than the first type 
of predictions, the latter ones have suffered in the past from larger uncertainties. 
In IP a QCD sum rule calculation of the strange quark mass in the scalar channel 
has been presented, which employed N^LO ((9(a^)) results for the correlator of two 
scalar strangeness-changing currents in perturbation theory (see also for a similar 
calculation) (for earlier calculations, see @, |^, ||, [l^, 

In the meantime the N'^LO {0{af)) correction to this correlator has become 
also known The perturbative contribution dominates the sum rule [P, 0, so 

one naturally expects the N'^LO correction to alter these results in a significant 
way. This expectation is also supported by a simple estimate of these corrections 
0] which shows that their omission is likely to constitute the main source of errors 
in the calculations of and 0. 

Since the N^LO correction is now known, we are in a position to present a 
reevaluation of the strange quark mass computation given in H]. 



2. Four- loop contributions to the scalar correlator 

The QCD sum rule used in this paper is based on the correlator of two scalar currents 
V'(Q',as,m„/i) = z / dxe'«"(0|TJ(s)jt(o)|0) (1) 



where J = daS'y°'u = i{ms — m„)sM, Q'^ = —q^ . It will be more convenient to 
work with the second derivative of 'ijj{Q'^), ijj"{Q'^) = d^/d(Q^)^, which satisfies a 
homogeneous renormalization-group equation fi-^ip" (Q"^) = 0. 

We will write ip"{Q) as ipp{Q)+ip'MpiQ)) "with ipp{Q) is the perturbative part and 
i'Npi.Q) contains the vacuum expectation values of the higher dimension operators. 
For the perturbative part one obtains the following result: 
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The terms of order in the 0(mq) part of 



have been extracted from the 
recent four-loop calculation of The terms of order in the O(m^) part of 
can be found in |2^ . 

The exact value of Oi agrees well with an estimate [0] of the same quantity^ 
based on the assumption of a continued geometric growth of the perturbative series 
for ip'^Q"^), which gave ai = 2660. 

We have neglected the light quark mass ruu, except in the overall factors. The 
renormalized parameters as and m^, m„ are taken at the scale /i. Their /i-dependence 
should cancel against that of the log /i factors in (|^) so that ip"{Q) is /^-independent. 

As for the nonperturbative contributions, we keep only the dimension-4 oper- 
ators. These are given, together with their renormalization-group properties and 
the values of the coefficient functions c, to next-to-leading order, in |p, ^ where the 
references to the original calculations can also be found. We quote here only the 
final result for ip'^p{Q): 
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Is and Ig are the vacuum expectation values of the two RG-invariant combinations 
of dimension 4, which are given for Uf = 3 and to the order we are working, by 



{ss)o + 



77r2 



TT 

as 



53\ 
24/ 



(4) 



^The corresponding constant in Q is called C31 — ai/6. 
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4. The sum rule 



To enhance the contribution of the low-lying states, one applies a Borel transform to 
the both sides of the dispersion relation used to define the sum rule 0. The effect 
is to transform the power-suppression of the states with a large invariant mass into 
an exponential one, controlled by the Borel parameter M^: 



1 1 



M6 TT 7o 

The Borel transform of the l.h.s. can be computed from 
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I) and is given by 
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The numerical constants entering these expressions have the values = —'-/e = 
-0.577 and ((3) = 1.202. 

At this point the usual procedure is to take advantage of the /i-independence of 
L[ip"{Q)] and the fact that the operation of Borel transformation does not act on 
/i) and choose /i = M |]2T|. This "renormalization-group improvement" effectively 
shifts the logs of M^//i^ into the renormalized parameters as{M) and ms{M). To 
the order we are working and for nj = 3, these are given by 
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We have denoted here L = \og{M'^ / Aq(jj^) and LL = log L. We have used in ([lOD 



the recently calculated exact values of the 4-loop beta function p2 
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and of the 4-loop mass anomalous dimension p5|, 
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The exact result (^) is close to an estimate based on the assumption of a geo- 
metrical growth for the coefficients of the 7^ anomalous dimension 

74(% = 3) = 4^^^^ 81.368, (13) 
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Table 1. Values of ml{M^){^^^^y at = 3 GeV^ (A^^ = 380 MeV) used for 
the discussion of the vahdity of the truncation approximation in the text. 



The usual approach followed in the numerical evaluation of the sum rule |^ 
has been to expand the Borel transforms (|^) and (§) in powers of l/L and truncate 
the resulting expressions to a given order in this parameter. There are certain errors 
inherent to this procedure which actually turn out to be important in practice. This 
can be seen at order 1/L^ by examining the structure of the Borel transform of the 
leading term in (|^ 

-■(-^) ^ IS^ (^) ^ - (^) ' — •) 

where we neglected the light quark mass. The first three coefficients Cj have the 
values ci = 4.8211, ca = 21.9765, C3 = 53.1421. 

We tabulated in Table 1 the values of the quantities m1{M^){as{M'^) /ny (with 
rfis = 1) for i = — 3 using two different approximations at a typical value of the 
Borel parameter = 3 GeV^. The first three columns show the values of these 
parameters computed by expanding in powers of the small parameter l/L up to the 
shown order. This approximation has been commonly used in the previous literature 
(for example the results in |^ have been obtained using a similar expansion to 
order 1/L^). 

From Table 1 one can see that truncating m^(as/7r)^ to order 1/L^ results in 
an error of the order of 100%, as the 1/L^ correction to its value is comparable 
with the 1/L^ term itself. The result of truncating to order 1/L^ is to overestimate 
the C(as) correction to the sum rule by a factor of 2. Neglecting this fact could 
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result in the paradoxical consequence that adding the C(a^) correction decreases 
the perturbative contribution Ti(M^) to the sum rule, although each separate term 
in the ct^ expansion is positive! 

The last three columns of Table 1 show the untruncated values for these pa- 
rameters (e.g. m^^\a[^^ means that the full i-loop expression for the running mass 
and the j-loop one for the running coupling have been used). One can see that 
these quantities are more stable when going from one order in perturbation theory 
to another than the truncated quantities. They can be therefore expected to give 
a closer estimate of the true size of each correction term and will be used in our 
numerical analysis below. 

Another difference from the treatment followed in [P, will be the use of the 
4-loop formulas for the running parameters ,|13) in all our expressions (|^ , and 
(plOl). We recall that the truncation approach (working to a finite order in l/L) 
employs running parameters at lower orders in the loop expansion in the power- 
supressed terms. The numerical differences between these two approaches are not 
significant. However, the former one is physically preferable as it uses as small ex- 
pansion parameter as{M) which is what is directly measured in practice (as opposed 
to l/L). 

The condensates entering the renormalization-group invariant condensates (^ 
and (^ are taken, as in 0, at the reference scale /xq = 1 GeV and have the values 
{uu)\^, = -(0.225)3 GeV=^ and (^G^) = 0.02 - 0.06 GeV^ The amount of SU(3)- 
breaking in the scalar condensate (ss) / (uu) will be varied between 0.7 and 1. The 
up quark mass has been taken as m„(l GeV) =5 MeV. 

The hadronic contribution to the sum rule is expressed below sq = 6 — 7 GeV^ 
in terms of the scalar form-factor d{s) in Ki^ decays. One has M 



^Im^(.) = ii-,\d{s)\\ 




(16) 



The values of s appearing in this relation are not accessible in Ki^ decays, but 
extend above the i^vr-production threshold Sth = {jtik + m.„Y. We will parametrize 
the scalar form-factor d{s) in this region by a sum of Breit-Wigner resonances cor- 
responding to the two bound states with the quantum numbers of the scalar current 
/s:*(1430) and K*{l^bQ) |TD|, |^ 



^ M(...)p(^^^--);+^™ (M|-.)^ + M|ri(.) ^^^^ 

1 1 7-12 
{Ml - {jUK + m^ff ^ {Mi - {niK + m^y^ 

The threshold value of the scalar form-factor d{sth) has been computed to order 
p'^ in chiral perturbation theory [jl[ with the result \d{sth)\'^ = 0.35 GeV^. The same 



7 



quantity has been recently extracted 0] from data on s-wave phase shifts for Ktt 
scattering [27| with a similar result \d{sth)\^ = 0.33 ± 0.02 GeV^. The values of 
the other parameters in (0) are 0] Mi = 1423 ± 10 MeV, Ti = 268 ± 25 MeV, 
Ma = 1945 ± 22 MeV, T2 = 201 ± 86 MeV, 7 = 0.5 ± 0.3. The energy-dependent 
widths Ti{s) are given by 
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For the region of large invariant mass of the hadronic states (s > sq), parton- 
hadronic duality can be expected to hold to a high degree of precision. This allows 
us to take the spectral density Im "^(s) to be equal to the imaginary part of the 
QCD expression (§). This is given by 
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The integration of this expression over t = {sq, 00) in (|^) is performed numerically 
keeping /i arbitrary, after which /i is set equal to M^ . At this point we have again the 
option of truncating or not the expression obtained after expanding ms(M), as{M) 
according to (|9|, p!0|) to a given power of 1/ ln(A2/M2). In accordance with the treat- 
ment of the theoretical side of the sum rule discussed above, we choose not to 
truncate the integral of the perturbative discontinuity either. 
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5. Results and discussion 



In Fig.l are presented plots of the invariant mass and of the running mass at 
the scale 1 GeV as a function of the Borel parameter for different values of the 
continuum threshold sq and the central value of the QCD scale = 380 MeV 

We extract our results from the region in M corresponding to the stability 



interval = 2 — 9 GeV^, obtaining in this way = 172 — 191 MeV respectively 
ms{l GeV)=191-213 MeV. The error arises mainly from the sq and dependence, 
the errors due to the condensates bein^ negligible, under 1-2%. 

The effect on rhg of changing between the limits 280-480 MeV is shown 

in Fig.2. The continuum threshold has been chosen such that optimal stability is 
obtained for each value of For A^^ = 280,380,480 MeV we find sq = 

5.0,6.0 and 6.9 GeV^. The corresponding values for the invariant mass are rn^ = 
231 - 232 MeV, 181-182 MeV and 140-147 MeV. This rather large spread of values is 
considerably reduced for the running mass at the scale 1 GeV, for which we obtain 
m,(l GeV) = 209-210 MeV, 201-202 MeV and 211-221 MeV. 

As one can see from Fig.2 the larger values of ^.^(l GeV) arises from including 
the large value of the QCD scale A^ = 480 MeV. If this curve is eliminated the 
following results are obtained: 

rhs = 181 - 232 MeV , m,(lGeV) = 201 - 210 MeV . (20) 



A similar observation has been made in p9[ in the context of the QCD sum rule for 



the p meson width, where even lower values for A^|^^ are advocated, of the order of 
220 MeV. We will adopt therefore in the following ( pOD as our result incorporating 
the theoretical errors arising from varying A^'^ = 280 - 380 MeV. 

These results are significantly higher than the 0{a1) results of [|, 0], so that an 
explanation for this difference is necessary. As mentioned already (see the discussion 
surrounding Eq.(16)) the approach used is this paper differs from that of [|], 0] in that 
the Borel transform is not expanded in powers of 1/ ln(A^/M^) but all orders in this 
parameter are kept. As an effect the leading term in the perturbative contribution 
to the sum rule (of (9(m^/M^)) is smaller than in @, |^, even after including the 
4-loop contribution, by about 8%. This results in an increase in the invariant mass 
by 4%, respectively 7-9 MeV. Another effect which pushes the result to the high side 
is the increase of the (9(mf/M'^) term, when adding the O^a"^) contribution. Since 
this term contributes with a negative sign to the theoretical side of the sum rule, it 
results also in a small increase of 1-2 MeV in the final result. 

For purposes of comparation we give also the results obtained if both sides of the 
sum rule had been truncated to order (in the leading terms of (9(m^/M^)). For 
A^l^^ = 380 MeV the best stability is obtained for sq = 5.7 GeV^ and the results 
for the strange quark mass are = 175 — 176 MeV, respectively 
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ms(l GeV)=192 — 193 MeV. Changing the continuum threshold sq by ±0.5 GeV^ 
about this value gives the broader range of values rhs = 167 — 184 MeV, ms{l 
GeV) = 183 - 202 MeV. Choosing = 280,480 MeV gives (for sq = 4.8,6.7 

GeV^) the results m, = 225 - 227, 135 - 139 MeV, respectively m,(l GeV)=203 - 
205, 204 MeV. These results are somewhat smaller than the ones obtained in the 
nontruncation approach (pO|) but still larger than the ones obtained in |^. The 
reason for this is that, as explained in Sec. 4, the leading term on the theoretical 
side of the sum rule decreases by about 6% when going from 1/L^ to This 
is partly compensated by a similar decrease in the contribution of the perturbative 
continuum when truncated to the same order in 1/L, such that the final result for 
rris is smaller than in the untruncated approach. 

Finally we should include the errors induced by the variation of the parameters 
(masses and widths) of the resonances and of the normalization factor d{sth)- The 
former give an additional error of about ±14 MeV on the value of rh^ and of ±15 
MeV on rhs{l GeV). The latter induces an error of about ±(11-12) MeV on both 
mass parameters. 

Adding all these errors in quadrature we obtain from (|20|) our final result 

m,(l GeV) = 205.5 ± 19.1 MeV . (21) 

This value lies on the high side of the existing QCD sum rule calculations of nis P, 
^, 0], coming closest to the recent result obtained to three-loop order in of 196. 7± 
29.1 MeV. The comparatively low results obtained in |^ were in good agreement 
with the lattice |TB[ results for the strange quark mass. With our new value the 
disagreement between the two is back in place. In |jl6[ the value rhs{2GeV) = 
128 ± 18 MeV was obtained which gives ms{lGeV) = 172 ± 24 MeV. Recently, 
a new lattice calculation has appeared |^ with lower results: in the quenched 
approximation rhs{2 GeV)=90±20 MeV, corresponding to rhs{l GeV)=121 ± 27 
MeV, and for rij = 2 an even lower value m<j(2 GeV) =70 ± 15 MeV, respectively 
ms{l GeV) =94 ± 20 MeV. Conceivable explanations for this discrepancy are a) 
significant systematic errors in the parametrization of the hadronic density and b) 
large contributions of direct instantons to the correlator of scalar currents |^ p^, |33[|. 
In our case, we consider b) to be little probable given the large scales = 2 — 9 
GeV^ at which our determination is performed (for an explicit calculation in the 
pseudoscalar current case see l3^). Thus further progress in improving the accuracy 



of the strange quark mass determination using the methods of the present paper 
can only come from a better knowledge of the hadronic density function. With the 
advent of a r-charm factory it should be possible to directly measure it in the future 
in semileptonic r decays. 
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Figure captions 



Fig.l Dependence of the invariant strange quark mass rhg on the Borel parameter 
and on the continuum threshold sq (the three lower curves). The upper three 
curves show the running mass at the scale 1 GeV for the same values of the param- 
eters. {Aqcd = 380 MeV) 

Fig. 2 Dependence of the results on the value of the QCD scale Aqcd- The contin- 
uous lines are the results for the running mass at the scale 1 GeV and the dotted 
lines show the invariant mass m^. 
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